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We present lattice results for the vector and scalar form factors of the semileptonic decays 
D — !■ Kivii and D K£v( in the physical range of values of squared four momentum transfer 
q^, obtained with Nf = 2 maximally twisted Wilson fermions simulated at three different lat- 
tice spacings (a ~ 0.102 fm, 0.086 fm, 0.068 fm) with pion masses as light as 270 MeV 
and m,tL > 4. The form factors are extracted using a double ratios strategy, which allows a 
good statistical accuracy and is independent of the vector current renormalization constant. The 
chiral/continuum extrapolation is performed through a simultaneous fit in the three variables 
(m;r,q'^,fl) using HMChPT formulae with additional ff{a^) terms that parametrically account 
for the lattice spacing dependence. Our results are in very good agreement with the experimental 
data in the full range for both D — !> niVi and D K£v(. At zero momentum transfer we obtain 
/o^'f (0) = 0.65(6)stat(6)syst and /^^^(O) 0.76(5)stat(5)syst, where the systematic error does 
not include the effects of quenching the strange and the charm quarks. Our findings are in good 
agreement with recent lattice calculations at A'/ = 2 + 1 . 
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1. Overview 

Weak decays of hadrons represent a very important source of direct information about the 
Cabibbo-Kobayashi-Maskawa (CKM) [||] matrix elements, which are fundamental parameters of 
the Standard Model (SM) flavour sector. The direct extraction of such quantities from the exper- 
imental decay rates requires theoretical inputs, namely the form factors and the decay constants 
which encode the non-perturbative QCD dynamics. Lattice QCD provides a way to compute with 
a good accuracy these quantities from first principles. 

In this contribution we present our determination of the form factors /^o^^(g^) and o^^(g^) 
obtained using the gauge configurations generated by the European Twisted Mass Collaboration 
(ETMC) adopting tree-level improved Symanzik gauge action and Nf = 2 twisted-mass lattice 
quark action, tuned at maximal twist [Q] to get automatic ^(a) improvement [||]. 

The present analysis updates and finalizes our previous study presented in ref. [Q], which was 
based mainly on simulations at a single value of the lattice spacing, a ~ 0.086 fm. Using the gauge 
ensembles ^ A2-A3 at j3 = 3.8 (a ~ 0.101 fm), B2-B4, and B(,-B-j at jS = 3.9 {a ~ 0.086 fm) and 
C2-C3 at j3 = 4.05 {a ~ 0.068 fm), we can now extrapolate safely to the continuum limit. The 
pion mass ranges from 500 MeV down to 270 MeV, and the size L of our lattices guarantees that 
m^L > 4, except in the case of the lightest pion for which we have trij^L > 3.7. For each pion 
mass and lattice spacing we have simulated several values of the (bare) strange and charm quarks 
mass to allow for a smooth, local interpolation to the physical values of and |^. We impose 
non-periodic boundary conditions on valence quarks |^, which enable us to inject arbitrary values 
of quark momenta in order to cover the full physical range. 

We extract the form factors using suitable ratios and double ratios of 2-point and 3-point 
functions (smoothly interpolated to the physical strange and charm quark masses) at fixed values of 
{m-n, q^,a). A combined analysis in these variables is then performed in order to reach the physical 
point. To this end, we fit our data with the predictions of SU(2) HMChPT [Q], which describe the 
mass and momentum dependencies originating from chiral loops in terms of a finite number of 
LECs, modified by the addition of ff{a^) terms parametrizing the lattice spacing dependence. 

Our preliminary results are in very good agreement with the experimental data in the full 
range for both D n£Vi and D KiVi. In particular at zero momentum transfer we obtain 
fD^^(0) = 0.65(6)stat(6)syst and /+^'^(0) = 0.76(5)stat(5)syst, where the systematic error does 
not include the effects of quenching the strange and the charm quarks. Our findings are in good 
agreement with recent lattice calculations at A^y = 2 + 1 [||], showing that the error due to the strange 
quark quenching is smaller than the present uncertainties. 

2. Formalism 

In the Standard Model, CKM matrix unitarity implies that flavour changing decays, at tree 
level, are due only to charged currents interactions. In the case of the semileptonic decay Hq^i 
Pqqi£V(, where HQqi is a heavy-light pseudoscalar meson and Pqqi is a light-light pseudoscalar 
meson (the subscripts Qq' and qq' represent the valence quark content') only the vector current 
contributes to the hadronic matrix element. The decay proceeds at quark level through Q — > q£v£, 

' 2 is a heavy quark in the sense that niQ 2> Aqcd, while q and q' are light quarks. 
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while the quark ^' is just a spectator. The standard Lorentz decomposition of such matrix element 
is done in terms of two form factors, f+{q^) and fo{q^), which encode the non-perturbative QCD 
dynamics: 

{P{k)\qr^Q\H{p))=f+{q^)[p^+k^-q^{ml-ml)/q'-]+fo{q^)q^{^^^ (2.1) 

where q^ = {k — p)^ is the four-momentum transfer, < < (m// — mp)^, and the two form 
factors obey the kinematical constraint /+(0) = /o(0). 

In the case of the semileptonic decay of a heavy meson, another convenient decomposition, in 
which the form factors are independent of the heavy meson mass in the static hmit, is given by 



{P{k)\qY^Q\Hip)) = y^{v^ + p^Jp{E)) (2.2) 

where v = ph/i^h is the //meson 4- velocity, p± = pp—Ev andE = v-pp = [m^+m^p — q^) / {Imn) 
is the P meson energy in the H meson rest frame. The two sets of form factors are related by 

/o(^2) = .J^/(ml-ml) [{mH-Ep)MEp) - plfp{Ep)] , (2.3) 
f+iq^) = [fv{Ep) + {mH-Ep)fp{Ep)] /mn . (2.4) 

The chiral and momentum behaviours of heavy meson form factors are described by the Heavy 
Meson Chiral Perturbation Theory (HMChPT). The formulae relevant for our analysis have been 
computed in continuum SU(2) HMChPT, at next-to-leading order in ChPT and at leading order 
in the HQET expansion in 1 / mh in the partially quenched and unquenched case with degenerate 
dynamical quarks [Q]. The explicit formulae for the form factors of the decay Hq^i — > Pg^i in 
our Nf = 2 lattice setup can be derived from the mentioned paper by computing the unitary limit 
= nisea- With the addition of a term accounting for discretization effects (starting at i^(a^) due 
to automatic ff{a) improvement in maximally tmLQCD), they read^ 

fp{m^'q',m,g,,E,a^) = + df^"'"'^'^'"'^ +C,{E)m],^, +C2{E) +€^0") , (2.5) 

/vKy, =£»o(l + 5/J''^''^''^''^+Di(£)m2,^,+D2(£)+D3a2) , (2.6) 

where the three contributions represent respectively the non-analytic {Sfp^) and the analytic terms 
(Ci,2, D\ 2), originating from HMChPT, and the discretization effects (C3, D3). The quantity 
A = m*fj — mii entering the pole factor in fp is the mass splitting between the vector H* and H. At 
leading order in the HQET expansion this splitting is zero, and in fact A is consistently neglected 
in all the loops. Nevertheless, it is customary to keep A in the tree level contribution to fp, since it 
correctly accounts for the position of the pole expected at q^ = mj^,. by vector-meson dominance. 

^We label as nixY the mass of a pseudoscalax meson with quark content xy. 
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For completeness, we collect here the explicit formulae for the non analytic terms: 

+ ng^ {-24ml^, - 9m2,^,m,, + 17m^^)/6£ , (2.7) 

+ 2h{m,,,,E) + (4,, -0^^^^/2, (2.8) 

which in the limit of degenerate valence quarks q' = q read 

(471/)^ 5/1;"^''^''''^ = -3{\ + 3g^)him,,)/4 + 4g^Mm,„E)-SKg^my3E, (2.9) 

{4nff5fj"'^'^'''"^ = (5 - 9/)/i(m,J/4 + 2/2K„£) , (2.10) 



where g is related to the coupling constant gDU'n- The functions /i, I2 and /i are defined in the MS 
scheme as in Ref. [Q] and the analytic terms depend on unknown LECs, functions of E. The range 
of applicability of the NLO chiral logs is generally expected to be limited to energies E <^ AchPT> 
which implies q^ close to ^^ax- However, using the so-called Hard Pion (Heavy Meson) ChPT 
it has been recently shown that the (HM)ChPT coefficients of the chiral logs are computable also 
at q-^ = 0, i.e. far from q^ ~ ^^^x- Thus, we will use these chiral predictions to parameterize the 
mass and momentum dependencies of lattice data in the whole q^ range. 

Following ref. [||], we are also working on a combined extrapolation based on the so-called 
z-expansion. The outcome of such analysis will be included in a forthcoming publication. 

3. Analysis 

The form factors are extracted from the lattice three-point correlation functions using a double 



ratios strategy [10, 11|, which allows a good statistical accuracy and is independent of the vector 



current renormalization constant, namely 

^^o(^max), (3.1) 
-^R\{q^), (3.2) 



plateau / 2 
C^/"{q,t) CP^{0,t)C""{0,t-T/2) pMeau 



Cf^(O,0 CPP{q,t)C""{q,t-T/2^ 

Cr"{Q,t) plateau , , 2\ o^ 

yR2{q), (3.3) 

where the ratios R'- are linear combinations of the form factors f+,o{q^)- 

The chiral/continuum extrapolation is performed through a combined {rriji^q^ ,a) fit using the 
modified HMChPT formulae ( ^ ) with a polynomial Ansatz for the unknown energy dependence 
of the analytic terms. We include in the fits the lattice data up to £" » 1 GeV and we allow terms 
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Figure 1: Pion mass (a), energy (b) and lattice spacing (c) dependencies of the combined fit for 
fpv^^ together with the statistical uncertainty band. The red points show the values of our fit, 
respectively, at the physical m-^, E corresponding to =0 and a = 0, while keeping the other two 
variables fixed. In (d) the results for f^^"" (q^) extrapolated to the physical point (together with 
the statistical error bands) are shown and compared with the CLEO-c experimental points /|7^. 



up to order E^ in the expansion of Ci{E) and Di{E). For the quantity g we adopt the most recent 
value g = 0.67(14) [p^, obtained hy Nf = 2 lattice simulations at fine lattice spacings, while 
the mass sphtting A is fixed to the PDG values, A(D tt) = 145 MeV and A(D K) = 
248 MeV. Finally, for the parameter / we choose the pion decay constant in the chiral limit /o = 
0.1215(l)(+o })MeV, as determined by the ETM collaboration [g]. 

The quality of our fit is illustrated in Figs. |l^|lc| for the D — > tt decay. Note that the discretiza- 
tion effects appear well described by terms of ff{a^), as expected. No appreciable dependence of 
the discretization effects on niji and E is observed. 
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— fo^q"^ t'lis work 

— f {q~) this work 

■ yq") CLEO-c D°->K 

A yq^) CLEO-c D*->K° 
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Figure 2: Results for the f {q ) form factors at the physical point versus q (a), compared 
to the CLEO-c measurements /|74|/. In (b) the same quantities divided by f^^^{0) = /(f^^CO) 
are compared to the experimental data of BABAR /173|7 and FOCUS /|7^ experiments. The bands 
represent the statistical error of the extrapolation. 



In Fig. and Figs. ^ ■ 2b we show the extiapolation of our lattice results to the physical value 
of mji and in the continuum limit in terms of f+{q^) and fo{q^). The bands represent the statistical 
uncertainty resulting from the analysis. 



3.1 Systematic errors 

The fitting procedure is subject to a number of sources of systematic uncertainties, which we 
are going to discuss in this section except for the effects of quenching the strange and the charm 
quarks. We present now a preliminary estimate of the systematic error at zero-momentum transfer, 
while our final estimate will be included in a forthcoming publication. 

Fitting function and energy range. We test the stability of our fits by adding/removing terms 
up to in the LEC's and/or possible NNLO corrections of order 0'{m%), as well as by includ- 
ing/excluding data with Eji > IGeV. Terms proportional to nij^E and E^ do not modify the fit of 
fp, which is essentially dominated by the pole factor, while they are necessary for describing the 
behaviour of /y. We estimate an overall uncertainty of 7% (5%) for D — > tt (D ^ K) decays. 

Discretization effects. Taking the difference between the result in the continuum limit and at 
our finest lattice spacing, we estimate that discretization errors are of the order of 5% (3%). 

Finite size effects. In the fitting procedure we include the lattice points with the lightest pion, 
mjc ^ 270MeV, having mj^L ^ 3.7. By excluding these data, we estimate that volume effects does 
not exceed 2% (2%). 

Value of the coupling constant g. We vary the value of the parameter g in the range [0.50,0.67] 
corresponding to available results in the Uterature. The uncertainty is below 3% (3%). 

Value of the mass splitting A. We also try to use the lattice determined value instead of the 
PDG one. This choice increases the statistical error as expected, since the vector meson mass is 
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poorly determined on the lattice, and it does not generate substantial variations in the central values. 
If treated as a free parameter, it is poorly determined. The estimated uncertainty is 1% (1%). 

In conclusion, at zero-momentum transfer we get the results 

f^^{0) = 0.65 (6)stat (6)syst, f^^iO) = 0.76 (5)stat (5)syst • (3.4) 

where the systematic error does not include the effects of quenching the strange and the charm 
quarks. Our findings are in good agreement with recent lattice calculations at A/y = 2 + 1 |^], 
showing that the error due to the strange quark quenching is smaller than the present uncertainties. 
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